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We show that Maxwell equations and Dirac equation (with zero mass term) have both sublu- 
minal and superluminal solutions in vacuum. We also discuss the possible fundamental physical 
consequences of our results. 
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I. INTRODUCTION 



According to Bosanac [Bo83| there is no formal proof based only on Maxwell equations that no electromagnetic 
wave packet can travel faster than the vacuum speed of light c (c = 1 in the natural units to be used here). Well, 
the main purpose of this paper is to show that Maxwell equations (and also the Dirac equation with zero mass) have 
superluminal solutions {v > 1) and also subluminal solutions {v < 1) in vacuum. 

This paper is organized as follows. In sect. 2 we introduce some mathematical tools that will be used. In sect. 3 we 
show how to construct the se called subluminal and superluminal solutions of the free Maxwell equations. In sect. 4 
the subluminal and superluminal solutions of the massless Dirac equation are discussed. Finally in sect. 5 we discuss 
some of the possible physical implications of these results. 
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II. MATHEMATICAL PRELIMINARIES 

In order to discuss these new solutions of Maxwell and Dirac equations in an unified way we briefly recall how these 
equations can be written in t he Clifford and Sp in-Clifford bundles over Minkowski spacetime. Details concerning 
these theories can be found in [|RS93yRS95|jRO90|l . 

Let A4 — {M, g, D) be Minkowski spacetime. (Af, g) is a four dimensional time oriented and space oriented 
Lorentzian manifold, with M ~ JR** and g g sec(r*M x T*M) being a Lorentzian metric of signature (1,3). T*M 
[TM] is the cotangent [tangent] bun dle. T*M ^ \J^(zmT*M and TM = U^^mT^M, and T^M ~ T*M ~ K^' , where 
R^'^ is the Minkowski vector space | |S W77| , |RR8"9| ] . D is the Levi-Civita connetion of g, i.e., Dg = 0, T{D) = 0. Also 
R{D) — Q,T and R being respectively the torsion and curvature tensors. Now, the Clifford bundle of differential forms 
CiiM) i s^the bundle of algebras C£{M) = VJ^^mCI{T*M), where Vcc e M,Ci{T*M) = Cii^s, the so called spacetime 
algebra | Lo93 IIS87]. Locally as a linear space over the real field R,C£{T*{M) is isomorphic to the Cartan algebra 
/\{T*{M) of the cotangent space and i\{T*M) = Y1=q f\^iT*M), where /\^{T*M) is the (^) dimensional space of k- 
forms. The Cartan bundle /\{M) — U^gM /\iT*M) can then be thought as "imbedded" in C£{M). In this way sections 
of C£{M) can be represented as a sum of inhomogeneous differential forms. Let {e^ = af?"} ^ secTAf, {fj, = 0, 1, 2, 3) 
be an orthonormal basis g(e^,ey) = 77^1. = diag(l, — 1, — 1, — 1) and let {j'^ = dx'^} € sec f\^{M) C secC£{M) be the 
dual basis. Then, the fundamental Clifford product (in what follows to be dennoted by juxtaposition of symbols) is 
generated by 7^7" + 7'^7^ = 27?^'' and if C e secC£(M) we have 



C = s + + ^^M-7^7" 



where 7^ 



dx^dx^dx^dx^ is the volume element and 



,bf,v,af,^p,p e secA (M) C secC£(M). For 



Ar e sec/\'^(Af) c secC{M),Bs G sec/\''{M) we define [ |Lo93|jHS87[ | Ar ■ = {ArBs)\r-s\ and Ar A B^ ^ (ABs)r+s, 
where ( )k is the component in /\^{M) of the Clifford field. 

Besides the vector bundle Cl{M) we need also to introduce another vector bundle C^spin (i,3)(-^) [Spin^(l,3) ~ 
SL(2,(r)] called the Spin-Clifford bundle. We can show that Ci'gpin^(i 3) (Af) ~ C£{M)/Tl, i.e., it is a quotient bundle. 
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This means that sections of C£spin^(i,3) (M) are some special equivalence classes of sections of the Clifford bundle, i.e, 
they are equivalence sections of non-homogeneous differential forms (see eqs.(]^,||) below). 

Now, as well known an electromagnetic field is represented by F e sec/\ (M) C secC^(M). How to represent 
the Dirac spinor fields in this formalism? We can show that even sections of C^sj,i„^(i 3) (M), called Dirac-Hestenes 

spinor fields, do the job. If we fix two orthonormal basis, S = {7'^} as before, and E = {7^ = R'j'^R — A(^7'^} with 



A[; g S0+(l,3) and R G Spin^(l,3), RR — RR — 1, and where is the reversion operator in C^i^a [Lo93,HS87], then 
the representations of an even section ip G secC-^spin^(i^3) (M) are the sections Vs and oiCi{AI) related by 



and 



1 



(1) 



(2) 



Note that tp^, has the correct number of degrees of freedom in order to represent a Dirac spinor field, which is not the 
case with the so called Dirac-Kahler spinor field. 

Let * be the Hodge star operator * : t\{M) ^^~^{M). Then we can show that if Ap e secl\^{M) C secC£(M) 
we have -kA = ^7^. Let d and 5 be respectively the diferential and Hodge codifferential operators acting on sections 



of A(^^)- If G sec A^(M) c secC£(M), then 5ujp ^ {-Y - 



d * ujp, with -k 



identity. 



The Dirac operator acting on sections of C£{M) is the invariant first order differential operator 



and we can show the very important result [Ma90|: 



+ d-^d~S. 



With these preliminaries we can write Maxwell and Dirac equations as follows [HeGG]: 

OF = 0, 



If m = we have the massless Dirac equation 



av^s = 0, 



(3) 



(4) 



(5) 
(6) 



(7) 



which is Weyl's one when ipYi is reduced to a Weyl spinor field. Note that in this formalism Maxwell equations 
condensed in a single equation! Also, the specification of ips depends on the frame S. 

If one wants to work in terms of the usual spinor formalism, we can translate our results by choosing, for example, 
the standard matrix representation of {7''}, and for -^j] given by eq.(|^) we have the following (standard) matrix 
representation: 



where 



Right multiplication by 



s - ibi2 bi3 - 1623 
-bi3 - ib23 s + ibi2 



-bo3 + ip ~boi + ibo2 
-box - ibo2 603 + ip 



(8) 



(9) 



gives the usual Dirac spinor. 

Before we present the sublumin al and superluminal solutions F< and we shall define precisely an inertial 
reference frame (irf) | SW77 , RR89 |. An irf / G secTM is a timelike vector field pointing into the future such that 
(/(J, /) = 1 and DI — 0. Each integral line of / is called an observer. A chart {x^) of the maximal atlas of M is 
called naturally adapted to / if / = d/dx°. Putting J = cq, we can find a — d/dx^ such that g(e^, C i,) = rj^^ and the 
coordinates (x^) are the usual Einstein-Lorentz ones and have a precise operational meaning |RT85|. x'^ is measured 
by "ideal clocks" at rest synchronized "a la Einstein" and a;' are measured with "ideal rulers" . 



2 



III. SUBLUMINAL AND SUPERLUMINAL SOLUTIONS OF THE MAXWELL EQUATIONS 



Let A S sec /\^{M) C secC^(M) be the vector potential. We fix the Lorentz gauge, i.e., d ■ A — —dA — such that 
F = dA ^ {d ~ 5) A = dA. We have the following: 

Theorem: Let H e sec/\^{M) C secC£{M) be the so called Hertz potential. If H satisfies the wave equation, i.e, 
d^U = ry^^a^a^H = -{dS + 6d)U = 0, and ii A = -5U, then F = dA satisfies the Maxwell equations dF = 0. 

The proof is trivial. Indeed A = -SU, then SA = -5^U = and F ^ dA = dA. Now dF = [d - 5) [d - 5)A = 
-{dS + Sd)A ^ Sd{Sn) = -S^dn = O since SdH ^ -dSU from d^U = 0. 

Now, since our main purpose here is to exhibit the existence of the new solutions we present only particular cases, 
leaving a complete study for another publication. To show the existence of a stationary solution Fq {dFg = 0) relative 
to a given inertial frame / with adapted coordinates (a;^) = {t,x^) introduce above, we choose the Hertz potential 

Ho(i, x) = (b{x) exp (75^)7172. (10) 

Since d^Ho = we get that 0(a;) satisfies the Helmholtz equation 

V0 + 17^0 = (11) 

The solutions of this equation are well known. Here we must comment the fact that the scalar wave equation has 



solutions which travels with speed less than c is known since a long time, being dis covere d by H. Bateman [Btl5| in 
1915 and rediscovered by several people in the last few years, in particular by Barut | BB92| ] [|. An elementary solution 
of eq.(pT|) with spherical symmetry is 

0(x)=C7 , r=^x^ + y^ + z^, (12) 

r 

where C is an arbitrary constant. From this result we can write Fq in polar coordinates (r, 9, tp) as 

C 

Fq = — [sin rit[a^r sin 9 sin ip — 13 cos 9 sin 9 cos (p)7°7"'^ 

— sin ^t{aVLr sin 9 cos + sin 9 cos 9 sin v)7°7^ 

+ sin m{f3 s\T? 9 - 2a)7"7^ + cos Vtt{(3 sw? 9 - 20)7^7^ 

+ cos rii (/3 sin 9 cos 9 sin p + aVLr sin 9 cos '^3)7"'^ 7'^ 

+ cos ilt{—l3 sin 9 cos 9 cos ip + ailr sin 9 sin 1^9)7^7^] (13) 

with a — fir cos fir — sinfir and f3 ^ 3a + fl^r^ sin fir. Note that Fq is regular at the origin and vanishes at infinity. 
Let us rewrite the above solution in terms of the old-fashioned vector algebra. We have that 



Eq = - sin fit, Bo = W cos fit, (14) 



where we defined 



W = C (^^- — P{x^ + y^) 2a \ ^^^^ 

\ ^3 ^3 ^5 ^3 / 

One can explicitly verify that divVF = 0, so that divii^o = divSp — 0, and that rotW^ + flW — 0, so that rot£^o + 
dBo/dt = and roti?o - dEo/dt = 0. 

We can show for a given F |Ma90 that 5° = ^F^^F is the 1-form representing the energy density and the Poynting 



vector. We have that Eq x Bq ~ 0, so that there is no propagation and the solution has vanishing angular momentum. 
The energy density is 

u= \[sm^9{fl^r^a^ +cos^ 9(3^) + {(3 sin^ 9- 2a)^] (16) 



^Barut found also subluminal solutions of Maxwell equations, with a procedure different from the one presented here, and his 
solutions are also different from the above. 
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Then J ud^x = oo. A finit e ener gy solution can be constructed by considering "wave packets" with a distribution of 
intrinsic frequencies f{fl) ||Btl5 |. These solutions will be discussed in another publication. It is also very important 



to see that Fq = Fq ■ Fq + Fq /\ Fq ^ 0, i.e, the field invariants are nonvanishing, differently of the usual solutions Fc 
of Maxwell equations that travel with constant speed c = 1 and for which F^ = 0. 

To obtain a solution moving with velocity < w < 1 relative to / it is necessary only to make a Lorentz boost 
in the x direction of the solution Fq. Another way to obtain a solution F'^ is to get a new solution $<(i,x) of the 
wave equation (9^<i>< = 0), like: 

^^(t,x)^C^-^^^eM^'i^t-kx)l (17) 
oj^-k^ = O^, (18) 

with 

^^^{jlix-vtr+y'+z^^^, 



7< 



1 dio 

VT^' '"~dk 



is a dispersion relation for a "particle" moving with velocity (here "group velocity" ) «< = dco/dk < I. From 
eq.(|l7|) we can obtain a new solution F^ moving with < v < I and which satisfies dF'^ = by taking F'^ = 9(9 •n<) 
with n< = $<7^7^. The fact that satisfies Maxwell equations follows from the above theorem, but it has been 
explicitly verified using REDUCE 3.5 g; unfortunately the explicit form of is very big to be given here. 
Now, to obtain a superluminal solution of Maxwell equations it is enough to observe that the function 

$>(t, x) ^ c '''''^^> exp[7^(cot - fcx)], (19) 

- fc2 = _f]2^ (20) 



with 



{jUx-vtf-y^-z^/^ 



1 du 

Vw^ - 1 "-k 



is a solution |BC93] of the wave equation c)^$> = 0, which travels with velocity v = duj/dk > 1. Writting n> = 
^>1^1^ we can obtain F> = d{d ■ n>) which satisfies 9F> — 0. This is then a superluminal electromagnetic 
configuration. The explicit form of Fy is again very big to be reproduced here, but again we verified it explicitly using 
REDUCE 3.5. Also F3 7^ which means that the field invariants are non null. We will discuss more the properties 
of this and some others extraordinary solutions of Maxwell equations in another paper. 

IV. SUBLUMINAL AND SUPERLUMINAL SOLUTIONS OF THE MASSLESS DIRAC EQUATION 

In order to find these kinds of solutions of the massless Dirac equation we shall make use of some ideas from 
supersymmetry. Since a Dirac-Hestenes spinor field tps is given by eq(y) we can define a generalized potential for tps ■ 
Indeed for each tps there exists A = A + j^B, A,Bg sec /\^ (A/) e secC£(M), such that 

= 9(^ + 7^5) (21) 



The use of REDUCE in Clifford Algebras is discussed, for example, in |Va95 Br87 
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The Dirac operator d plays here a role analogous to that of supersymmetry operator | Gr94[ . In fact, Clifford algebras 
are Za-graded algebras, and C£+C£+ C Ci+ , Cl^Ct^ c C^", C£-Cr c C€+, whereCF [Ci'] denotes the set of 
elements of CI with even [odd] grade. Since the Dirac operator has vector properties, its action transforms fields of 
even grade into fields of odd grade, and vice-versa. Representing a spinor field by means of nonhomogeneous forms 
of even degree is therefore equivalent to find a bosonic representation of a fcrmionic field. 

The quantity A can be interpreted as a kind of potential for the massless Dirac field. Indeed, from eq.(0) with 
eq.(|2l|) it follows that 



d^A = 0, 



(22) 



or that 



d'^A = 0, d'^B = 0. 



(23) 



A simple subluminal solution at rest relative to the inertial frame / = eo in the coordinates {x^) naturally adapted 
to / is 



with 0(5;) given by eq.(O). We have 



C 



^o(t, x) = 7°(/'(f) exp {-t^nt) 



(sin rir cos VLt'y^ — sin Vlr sin Vlt'-^^^f^^'^) 



(24) 



(25) 



Then 



C 

= — [-Jlr^sinilrsinilt 

-l-7°7^ Ace cos fit + 7°7^Ay cos fit 
-|-7"7^Azcos fit - ^^^"^Xz SYR fit 
+7^7^ Ay sin fit - 7^7^ Aa; sin fit 
-|-7°7^7^7^rir^ sin fir cos fit] , 



(26) 



where A = fir cos fir — sin fir. 

The above solution in the usual formalism reads 



/ 



i sin fit 



\z fl . ^ 
— r- + i — sm s2r 



J., 

i sin fit 



r 

- iy 



A 



:xz fi 

cos fit [ — + i — sin fir 

r 



V 



r 

cos fit 



X 



ly 



A 



(27) 



and one can explicitly verify that indeed dip'~' = 0. 

Other subluminal solutions can be obtainned by appropriated boosts. An explicit superluminal solution ijj^ can be 
obtained by writting 



A>{t,x)^j"<^>>{t,x) 



(28) 



with $> given by eq.(|l9|) and = [d^y{t, x)]'-f^. Agai n the explicit form of '5^ is very big and will be not written 
here, but it has been verified using REDUCE. Barut [ Ba94 has found subluminal solutions of d'i> = 0, where 
is the usual Dirac spinor field, with a different method from the one used here; our one is much simple since it is 
representation free and uses elegant tools from supersymmetry. 



5 



V. CONCLUSIONS 



We want to discuss three possible implications of the results we have shown. 

(i) If the superluminal solutions of at least Maxwell equations are realized in Nature we can have a breakdown of 
Lorentz invariance. Indeed, suppose / — d/dt is the fundamental reference frame and /' = — V^)d/dt — 

y/Vl — V^d/dx is the laboratory frame (an inertial frame). Suppose that Fy is a superluminal solution of Maxwell 
equations, i.e., i9F> = (w'^ — k'^ = — fi^), travelling forward in time according to /. Then, the validity of active 
Lorentz invariance implies that there exists R € Spin_|_(l,3) such that F!^ = RFyR satisfies dFi, — with going 
backward in time (and carrying negative energy) according to /. This solution can be interpreted as an "anti-field" 
coming from the past, and is a good solution. However, the physical equivalence of all inertial reference frames implies 
that according to /' there exists solutions Fy of Maxwell equations travelling forward in time and carrying positive 
energy E' — u)' {h = 1) according to /' but travelling backward in time (and carrying negative energy) according to 
/. The field F^ can be absorved, e.g., by a detector in periodic motion in / (it is enough that at the time of absortion 
the de tector has relative to / the velocity V of the /' frame). This generates as is well known a causal paradox 
(Tolman-Regge paradox). The possible solution is to say that / and /' are not physically equivalent. We then 



have the following: /' cannot send to some observers (integral lines) of the / reference frame a superluminal signal 

When w' ^ , . ^ 

We 



such that oj' < (V/Vl — V'^)il. When ui' = {V/\/l — V'^)n the superluminal generator of /' stops working for k' in 
some spacetime directions, and an observer in / can calcu late his abso lute velocity which is ^ = lo' j \fu/^- 



must also call the reader's attention that recently Nimtz | IIN94 , EN93 | transmited Mozart's symphony # 40 at 4.7 c 
through a retangular wave guide, that as i s now well known acts like a potential barrier for light. Important 

related results have also been obtained in |St93|. We can show easily that under Nimtz experimental conditions the 
solution of Maxwell equations in the guide gives a dispersion relation like eq(po|), i.e, corresponding to superluminal 
propagation |Re86|. We shall discuss this issue in details elsewhere. 



(ii) The existence of the subluminal solutions are very important for the following reason: Recently | VR93 , VR9f: 1 
we proved that dF = for F'^ 7^ is equivalent under certain conditions to a Dirac-Hestenes equation dipsj^j'^ + 
n^ip'sl" — 0, where F — V's7^7^'0s- This means that eventually particles are sp ecial stationary electromagnetic waves 
and a de Broglie interpretation of quantum mechanics seems possible [ RV93 1 . We will discuss this issue in details 
elsewhere. 

(iii) Finally the existence of subluminal and superluminal solutions for dtpY: = (which reduces to Weyl equation for 
'i/'E a Weyl spinor) may be important to solve some of the mysteri es associated with neutrinos. Indeed if neutrinos 
can be produced in the subluminal and superluminal modes - see [ pt95| , |Gn86|] for some experimental evidences for 
superluminal neutrinos - then they can eventually escape detection on earth after leaving the sun. Moreover, for 
neutrinos in a subluminal mode it would be possible to define a kind of "effective mass" . Recently some cosmological 
evidences that neutrinos may have a nonvanishing mass have been discussed [Pr95|. One such "effective mass" could 



be responsible for those cosmological evidences, and in such a way that we can still have a left-handed neutrino since 
it would satisfies the Weyl equation. We are going to study this proposal in a forthcoming paper. 
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Note added in proof. After we finished this paper we have been informed by Professor Ziolkowski that he and 
collaborators found also superluminal solu tions of the scalar wave equation and also of Maxwell equations and even 
Klein-Gordon equation [DZ92,DZ93, Zi93|. Also D r. Lu and collaborators found a very interesting "superluminal" 
solution of the scalar wave equation |LG92 .LG92a| and Lu even realized an approximation for his solution (the so 
called X waves) as a nondispersive pressure wave in water which travels with velocity 1.002 c, where c here is the 
velocity of sound in water! 

The solutions found by Lu can be used to construct Hertz potentials for the Maxwell equations and then to generate 
superluminal electromagnetic field configurations. We believe that it is in principle possible to build such fields with 
appropriate devices. This is also the opinion of Dr. Lu. Also some of the Ziolkowski solutions, according to his 
opinion, may be realized in the physical world. We shall discuss these points in another opportunity. 
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